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Global and local relaxation of a spin-chain under exact Schrodinger and 

master-equation dynamics 
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We solve the Schrodinger equation for an interacting spin-chain locally coupled to a quantum 
environment with a specific degeneracy structure. The reduced dynamics of the whole spin-chain 
as well as of single spins is analyzed. We show, that the total spin-chain relaxes to a thermal 
equilibrium state independently of the internal interaction strength. In contrast, the asymptotic 
states of each individual spin are thermal for weak but non-thermal for stronger spin-spin coupling. 
The transition between both scenarios is found for couplings of the order of 0.1 x AE, with AE 
denoting the Zeeman-splitting. We compare these results with a master equation treatment; when 
time averaged, both approaches lead to the same asymptotic state and finally with analytical results. 
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I. INTRODUCTION 

Various attempts have been made to account for 
thermodynamical behavior of quantum systems jj, 0- 
Especially the relaxation into an equilibrium state in 
a pure quantum world does not seem to be feasible 
since Schrodinger dynamics is reversible like the classi- 
cal Hamiltonian dynamics. By introducing irreversibility 
into quantum mechanics one thus has to face all the old 
difficulties. 

One quite successful way to introduce relaxation be- 
havior into quantum mechanical models is to consider 
open systems modeled by a quantum master equation 
0, or the Lindblad formalism 0, IE 13 ■ I n these ap- 
proaches the influence of the environmental system enters 
the Liouville-von-Neumann equation for the considered 
system via incoherent damping terms. To deduce such 
a closed evolution equation several approximations are 
necessary, e.g., the Born-Markov assumption. 

In the context of quantum thermodynamics has re- 
cently been found irreversible behavior in classes of very 
small bipartite quantum systems described by a pure 
Schrodinger evolution only 0, 0- This approach does 
not need those specific assumptions about the environ- 
ment. Instead, system and environment are treated as a 
whole. Even in a small bipartite systems, consisting of a 
two level system ("gas system") coupled to an environ- 
ment ( "container" ) of no more than some hundred levels 
thermodynamical behavior is generic - relaxation occurs 
to a theoretically predicted equilibrium state (see 01). 

In the above mentioned scenarios the gas system has 
been very small (from two to five levels) and coupled to 
an environment without any structure or selectivity. In a 
more complex situation the system under consideration 



(gas system) could be constructed from several identical 
subsystems, e.g. the system could be a spin chain. For 
such bipartite systems with increased internal complexity 
the effect of the coupling topology is not yet completely 
understood: If we couple a chain of identical subsystems 
at one edge only to a quantum environment as before, the 
question arises whether the whole system will still relax 
into a thermal equilibrium state. Furthermore, we are in- 
terested in whether and when the individual subsystems 
are also in a thermal equilibrium state. 

According to recent investi gati ons on local tempera- 
ture of modular systems [ll], H^h we expect the same 
global as well as local temperature in the system for 
weakly coupled chains. In cases of a stronger coupling 
this may no longer be the case (as shown by an open 
system approach 0). 

Those chain systems coupled to a quantum environ- 
ment will be treated as a closed system subject to 
Schrodinger dynamics and will be compared to the mod- 
eling via a quantum master equation. In particular, 
we will analyze chains with different coupling types and 
strengths with respect to their local as well as global 
thermal or non-thermal properties. 



II. THEORETICAL BACKGROUND 

A. The Considered System 

We consider chains of three identical quantum systems 
described by the Hamiltonian 

3 .2 

4 = ^#ioc(m) + 7 ^ s (m,m + i)- (i) 
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Here the first sum contains the local Hamilton operators 
H\ oc (/j) of site. In our case the local Hamiltonian of site 



2 



\x reads 



with n = 1,2,3 (2) 



where a z denotes the Pauli spin operator. Here and in 
the following all energies are taken in units of the Zceman 
splitting. The second sum refers to the next neighbor 
couplings between the subsystems normalized by 



a way that it couples only to one boundary of the system 
(see Fig. ^| , say to the third spin of the chain. Further- 
more, the interaction should allow for energy exchange 
between system and environment. 

In the following we will consider both, a complete so- 
lution of the Schrodinger equation for system and envi- 
ronment according to fftot as well as a quantum master 
equation, simulating a corresponding bath coupling. 
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Tr{/ S 2 } , 



(3) 



C. Hilbert space average approach 



where n is the Hilbert space dimension of the chain. 
Thus, it is possible to control the internal coupling 
strength by the single parameter A (see l[T]l) only, irre- 
spective of the type of coupling. In the following we 
compare two different interactions - a random coupling 

3 3 

%(p,H + l) =^X)ptf0-<(/i)<8>Oj-O* + l), (4) 
i=i j=i 

where the pij are normally distributed random numbers 
in the interval [—1, 1], and a Heisenberg interaction, 
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B. The Environment 
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FIG. 1: Coupling topology of the spin chain to the quantum 
environment. 

The model system described in the last section, is now 
coupled to an adequate environment: 



-fftot, — ffs 



ffe + y — ^si 



(6) 



The first term is the Hamiltonian of the chain, H e refers 
to the Hamiltonian of the environment and I se is the 
interaction between system and environment. After nor- 



malizing J se via ige = Tr{/ S 2 c } (compare Q) the ex- 
ternal coupling strength is controlled by the parameter 
k only. This system environment coupling is taken to be 
small in all cases as a precondition to allow for thermo- 
dynamic behavior. The interaction will be chosen in such 



Let us start with the complete model under 
Schrodinger dynamics. At first sight it may not seem 
clear how a totally time reversible equation like the 
Schrodinger equation could produce something like ir- 
reversible behavior. To clarify this point we introduce 
some aspects of quantum thermodynamics. A complete 
derivation of all aspects of this theory is beyond the scope 
of this article. Therefore we only sum up some central 
aspects of this approach and refer the interested reader 

to nm. 

According to quantum thermodynamics (see 9]), the 
system proper relaxes to a Gibbsian state whenever the 
density of states of the environment is an exponential 
function of energy; the latter is typical for many body 
systems. We model such a quantum environment here by 
a system (Fig. ^ the width of the energy levels should 
indicate the degeneracy) of the eight energy levels Ef 
with degeneracies given by 



N e {E° 



A£2 B * 



(7) 



This may seem rather artificial, but the idea is as follows: 
the environment may possess a continuous spectrum, but 
due to the weak coupling between s and e the system 
only couples to the resonant levels. Thus we can neglect 
all the other levels. This environment is then expected 
to induce on the spin-chain a canonical state with the 
reciprocal temperature /?, 



k B T 



d_ 

dE 



In N C (E) 



(8) 



where k-Q is the Boltzmann-constant. With the special 
degeneracy structure iQ this leads to (3 — In 2. 

For a numerical test and to avoid any bias we pick for 
I sc a hermitian random matrix from an ensemble with 
the distribution 0] 



P(Ise) = 



47T 



(9) 



Theoretical Predictions: In contrast to the master 
equation approach, the dynamics of the closed system 
(spin chain and environment) is purely Schrodinger-type. 
The von Neumann entropy of the completely closed sys- 
tem thus remains constant. But by splitting up the whole 
system into two parts, a small system taken to be the 



3 



spin-chain and a large one, the environment, the von 
Neumann entropy of the parts can change in time. As a 
matter of fact, the small system shows a thermodynam- 
ical behavior if two restrictions are met: 

a) The coupling between the small system and the 
environment should be small, i.e. the energy con- 
tained in the interaction has to be much smaller 
than the local energies, 



/so > < < H, 



H, 



which guarantees that the spectrum of the environ- 
ment is not disturbed too much. 

b) The Hilbert space of the environment should be 
very large compared to the Hilbert space of the 
small system. 

Conversely, all systems meeting these conditions may be 
called thermodynamical. 

The concrete form of the interaction defines an accessi- 
ble region within the whole Hilbert space, selected by the 
initial state. Using a model with full energy exchange be- 
tween system and environment (canonical situation) only 
one supplementary condition remains active, the overall 
energy conservation defining the accessible region. By 
a topological investigation of the Hilbert space (see |(|), 
the details of which are beyond the scope of this article, 
it is possible to show that the state of the complete sys- 
tem will enter a very large region (the dominant region) 
within the accessible region, for which the system under 
consideration is in a state with approximately maximum 
von Neumann entropy. The respective energy distribu- 
tion of the small system in this dominant region, i.e. the 
probability to find the small system in an energy eigen- 
state E® (no degeneracy) for an initial state with a sharp 
energy is then given by 



W d {E s i 



N tnt 



(11) 



(for a complete derivation of (| 1 1 ft see 0,^3)- * s the 
i-th energy level of the spin system, Ef the correspond- 
ing energy level in the environment with the degeneracy 
Ni(Ef) and iV tot = Y%=i N i( E t) is the total number of 
levels in the environment. 



The entire dynamics of the system coupled to the bath 
is given by the Liouville-von-Neumann equation, from 
which the Nakajima-Zwanzig equation can be derived 
Assuming that the bath is in a thermal state, 



Pc 



(12) 



an expansion of the latter up to second order in the sys- 
tem bath coupling reads: 



(10) dp s 



dt 



H s ,p s 



- I dt' Tr e 
o 



I sc (t' -t),p s 



where p s the reduced density matrix of the spin-system 
and I se (t), the system bath interaction in the interaction 
picture reads, 



I sc (t) = e 



-i(H a +H c )t 



(14) 



the time dependence of which can be computed in the 
eigenbasis of H s and H c . Applying, as usual, the Markov 
approximation, the integral of (| 1 31) can be computed and 
the following form of the damping rates is found: 



dp s _ 
dt 

where 



H s ,p s 



-(Ar a p s )-(r a p s A)+(p s r b A)-(Ap s r b ) 

(15) 



A = i(i) ®i(2) ®& x (3). 



(16) 



As can be seen the environment has been locally coupled 
to spin 3 only. 

Denoting the energy eigenvalues and eigenvectors of 
the spin-system by E l s and \i), respectively 



U) = Ef \i) 



(17) 



and defining = h(E? — Ej), the transition matri- 
ces r a / b have the following matrix representation in the 
eigenbasis of H s : 



,^-/3c-l' Ai J' 



<*|r Q \ j) = K 



A.. 



for ujij > (18) 



for Uij < (19) 



D. Master equation approach 

Master equation approaches describe the dynamics of 
open systems. They have been widely applied to describe 
system bath models, in particular in quantum optics 0- 
Their derivation is standard and can be found in several 
textbooks. However, since our system, the spin chain, has 
an internal structure and since only one of the boundary 
spins directly couples to the bath, special care must be 
taken: 



(i\T b \j) =k ^"l^ Aj, for u ij >0 (20) 



(*l rb |j) - v^-i - 4 ^ for ^' <0 (21) 

Note that by virtue of (|14|l the real damping rates of l|18(l 
- Ij21|) only appear if the transition matrices T a ^ b are rep- 
resented in the eigenbasis of H s . Therefore, simply writ- 
ing down Lindblad damping terms |l6j for spin 3, the one 
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which directly couples to the bath, would in general, lead 
to wrong results [17]]. However, if the coupling between 
the spins in the chain is small enough, additional approx- 
imations can be made and one will obtain the Lindblad 
type damping rates for spin 3 |18|. 

It should also be mentioned, that the thermal state 



(22) 



is the stationary and therefore asymptotic solution of 

03- ' 



III. NUMERICAL RESULTS 

A. Spectral-Temperature 

We will use as a measure to characterize the asymptotic 
state of the spin system and the single spins a " spectral- 
temperature" defined as 



1 



:=- 1 



W + Wj 



M 



M 

E 

71=1 



w n + w n 



lnVF»-lnW» 
E„ — E n -i 



(23) 



where W n is the probability to find the system at the n-th 
energy level E n , M is the highest and the lowest energy 
level. The main idea here is to assign a Boltzmann factor 
to each pair of neighboring energy levels. The spectral- 
temperature then simply is the average over all these 
factors weighed by the corresponding occupation prob- 
abilities. This spectral-temperature also exists in non- 
equilibrium situations but coincides with the standard 
definition of temperature (e. g. [Tflj l'l only for a canonical 
state. 

For the calculation of the spectral-temperatures ac- 
cording to (|23|) we insert for W n the long time average 
(Ps n (t)) t 01 the respective occupation probability. 



B. Schrodinger dynamics 

To demonstrate the relaxation behavior of the spin- 
chain under Schrodinger dynamics we have diagonalized 
the total Hamiltonian H tot (see (0) and solved the exact 
Schrodinger equation 



di 



l&ot(*)> 



:#tot |Vtot(*0 = 0)) 



(24) 



The initial state IV'tot(^o)) bas been taken as a product 
state of the spin chain and the environment. From the 
density matrix ptot(t) — \i>(t)) we trace out the 

environment and transform the resulting reduced density 
matrix into the eigenbasis of the spin-system. The time 
evolution of this reduced density matrix p s (t) is plotted 



in Fig. |21 for an internal random coupling and in Fig. [3| 
for a Heisenberg coupling. (We restrict our analysis to 
A > 0, the antiferromagnetic case. Corresponding results 
can be obtained for A < 0.) 




FIG. 2: Relaxation of the spin chain s into equilibrium un- 
der Schrodinger dynamics for A = 0.4 and k = 0.001 with 
a random interaction. The horizontal lines are the expected 
probabilities from 11 U to find the system at the correspond- 
ing energies. These are in good accordance with the time- 
averaged values of p B (t). A and k are energies taken in units 
of the Zeeman splitting. (Note that W d {E%) = W d {E%) ltTB.1 
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FIG. 3: Relaxation into equilibrium under Schrodinger dy- 
namics as of Fig. [21 but with a antiferromagnetic Heisenberg 
interaction. 

Both evolutions have been obtained with the coupling 
parameters n — 0.001 and A = 0.4. The horizontal lines 
are the probabilities as expected from Qllfl. As can be 
seen, for both internal couplings, the spin-system relaxes 
into a state which is in accordance with {TJ . The proba- 
bility fluctuations can be interpreted as a finite size effect 
of the environment. 

Global and local Temperatures: Now that we have 
demonstrated the relaxation of the spin-chain into a ther- 
mal equilibrium state, we want to analyze the relaxation 
behavior of the spin-chain depending on the internal in- 
teraction strength A. The approach of III CI is indepen- 
dent of the internal structure of the small system. Thus 
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the relaxation behavior should be independent of A and 
the temperature of the spin-chain should be identical to 
that induced by the environment T° = ^ (see (0)). To 
specify the temperature of the spin-chain T s we use the 
"spectral-temperature" defined in l|23|l . We compare the 
spectral-temperature of the total spin-chain to those of 
each individual spin. For the local spectral-temperature 
of a single spin, l|23[> reduces to 



T^IOC 

A* 



fr-rloc 



F^ - F^ 



ln(Wf) - ln(W ") : 



(25) 



where W/ 1 is the long time average of the probability to 
find the /x-th spin at the energy . These quantities 
are plotted as a function of A for a random interaction 
in Fig. 0] and for an antiferromagnetic Heisenberg inter- 
action in Fig. [5] 
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FIG. 4: Global and local temperatures as a function of A for 
a random interaction. 



more and more with increasing A although the whole 
spin-system continues to reach a canonical state 12]. 

Especially for a Heisenberg coupled spin chain one 
can verify analytically that for a canonical state with 
some respective temperature, the local temperatures de- 
viate more and more with increasing A (see Sec. IIVII . 
Note, that for a ferromagnetic Heisenberg coupling the 
spectral-temperatures of each spin show a different be- 
havior: they decrease with increasing A. In any case, 
they deviate from T c . 

To verify that the state of the total spin-chain p s is 
indeed canonical we test the off-diagonal elements. All 
absolute values are smaller than 10~ 4 . Therefore we ar- 
gue that the state of the total spin-chain is indeed a 
canonical one for all practical purposes and its spectral- 
temperature can be identified as the thermodynamic tem- 
perature. 

The deviation of the local spectral-temperatures from 
the temperature of the whole spin-chain can also be un- 
derstood by analyzing the correlation C between the 
spins 



C = Tr 



{[<Pi(*)>t® (Pa(*)> t ® (Ps(t)) t - (Ps(t)) t } 2 }. 



(26) 

The brackets () t denote the time average of the density 
matrices of the single spins (p M ) and the spin system (p a ). 
Due to increasing A the spins are more and more corre- 
lated. This correlation causes an increase of the local en- 
tropy which leads to an increased spectral-temperature 
of each individual spin. Fig. shows C as a function of 
A for a random interaction and Fig. [7|for a antiferromag- 
netic Heisenberg interaction. 
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FIG. 5: Global and local temperatures as a function of A for 
a antiferromagnetic Heisenberg interaction. 

The spectral-temperature of the spin chain T s ap- 
proaches the temperature T c imposed by the environ- 
ment, irrespective of the internal coupling strength A. 
However, the spectral-temperatures of each single spin 
is found to increase with increasing A. The reason for 
this behavior is that the local spin-system is disturbed 
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FIG. 6: Correlation C 1261 1 as a function of A for a random 
interaction. 

As can be seen from both figures, the correlations C in- 
crease with increasing A. The ferromagnetic Heisenberg 
coupling also show increasing correlations in the consid- 
ered ranges of A and T c . 



6 



0.12 



0.07 - 




0.02 1 1 1 1 1 1 1 1 1 

0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

A [AB S ] 

FIG. 7: Correlation C 1261 as a function of A for a Heisenberg 
interaction. 
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FIG. 9: Time evolution of p a under a master equation with 
antiferromagnetic Heisenberg interaction, (compare Fig. |3J. 



C. Master equation 

We have solved the master equation of Sec. Ill Dl with 
the same values for the parameter A = 0.4 and k = 0.001 
as for the Schrodinger dynamics. The corresponding re- 
laxation into equilibrium is shown in Fig. for a random 
coupling and in Fig. [5] for the antiferromagnetic Heisen- 
berg coupling. 
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FIG. 8: Time evolution of p s under a master equation with 
random interaction. The equilibrium reached is the same as 
that under Schrodinger dynamics (Fig. and the one pre- 
dicted by lllll (horizontal Lines). 



calculated the spectral-temperature H24JI for the whole 
chain and for each individual spin 1)25(1. The results are 
plotted in Fig. for a random interaction and in Fig.llll 
for a antiferromagnetic Heisenberg interaction. 
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FIG. 10: Spectral-temperatures under a master equation 
with random interaction as function of A. The solid line is 
the spectral-temperature T B of the total spin chain. The 
dashed line is the local temperature Tl° c of spin 1, the nar- 
rower dashed line T^° c the one of spin 2 and the dotted line 
Tj oc the one of spin 3. 



The horizontal lines denote, again, the equilibrium 
state, which the spin-system should reach according to 
pi[l. The state obtained via the master equation with 
only one spin coupled directly to the bath approaches 
the same equilibrium state as the Schrodinger dynamics. 
Also the off-diagonal elements are damped away. Asymp- 
totically both methods show, under the analyzed condi- 
tions, the same behavior whereas the short-time behavior 
is difficult to compare. 

Global and local Temperatures: As for the Schrodinger 
dynamics we now study the global and local states of the 
spin chain in dependence of the parameter A. We have 



In both cases the temperature of the whole spin-chain 
T s reaches the same value as the temperature of the 
environment T c = 1.44. For weak spin-spin coupling 
(A <C 1) the local spectral-temperature of each spin is 
approximately the same as the one of the whole chain. 
As for the Schrodinger dynamics with increasing A the 
local spectral-temperature of each spin rises. Especially 
the local temperature of the spin in the middle of the 
chain, spin 2, increases more rapidly than those at the 
boundaries. 
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FIG. 11: Spectral-temperatures under a master equation 
with antiferromagnetic Heisenberg interaction in dependence 
of A (compare Fig. I1UI . 



IV. ANALYTICAL RESULTS 

Finally we want to compare the numerical results 
shown before with the analytic solutions for the Heisen- 
berg coupled spin-chain. We start from a density-matrix 
ps{f3,\) for the whole spin-chain in a canonical state. 
Transforming ps((3, A) into the product basis, one can get 
the reduced density matrices of the single spins A). 
We have calculated the local spectral-temperature for 
each spin using (|23[) in dependence of (3 and A. The result 
is shown in Fig. 1121 (with the same global inverse temper- 
ature (3 = In 2 as for the numerical calculations) and is 
in good agreement with our numerical results (compare 
with Fig. |3]and Fig. ^] respectively). 




-0.4 -0.2 0.2 0.4 

FIG. 12: Analytical results for local spectral-temperatures 
with the global inverse temperature (5 — In 2 as function of 
A for ferromagnetic (A < 0) and antiferromagnetic (A > 0) 
Heisenberg coupling. 

To analyze the analytical solution for the correlation- 
measure C(/3, A) we have used again (|2"5)l with the an- 
alytical density matrices ps(/3,\) and A). Fig. H31 
shows C((3, A) for the global inverse temperature (3 = In 2. 
Again the analytical result is in very good accordance 
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FIG. 13: Correlations C as function of A (/3 = In 2) for ferro- 
magnetic (A < 0) and antiferromagnetic (A > 0) Heisenberg 
coupling. 

with our numerical results (compare with Fig. UJ. An 
interesting point is, that the correlations C also increase 
in the ferromagnetic case although it is known that there 
exists no nearest neighbor entanglement poj ]. 

Fig. ^] shows C(/3, A) as function of (3 and A. C in- 
creases with increasing (3 and A in the antiferromagnetic 
case as one would expect because of the entanglement in 
the spin-chain for low temperatures. For the ferromag- 
netic coupling one can see a bump depending on f3 and A. 
For low temperatures (large f3) C decreases and vanish 
for T = 0. 

Presently it is not quite clear why this bump occurs in 
the ferromagnetic case for intermediate temperatures. It 
can be shown, that the crest of the ferromagnetic bump 
is a function of (3 and A. The behavior for high temper- 
atures is the same for the ferromagnetic as well as the 
antiferromagnetic case: C is vanishing because the local 
states and the global state will be a totally mixed one. 



C 




A[AE S ] 



FIG. 14: Correlations C as function of /3 and A for ferro- 
magnetic (A < 0) and antiferromagnetic (A > 0) Heisenberg 
coupling. 

Note that C is not an entanglement measure. C checks 
whether ps(j3,X) is factorisable by the reduced density 
matrices Pfj,((3, A). Therefore C > in the ferromagnetic 
case for T > indicates on one hand that the product 
basis differs from the eigenbasis of the spin-chain 
(which is a pure quantum mechanical effect) and oth- 
erwise whether a local thermodynamical description is 
possible. 
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V. SUMMARY AND CONCLUSION 

The main motivation of this paper has been to an- 
alyze the thermodynamic behavior of one part ("the 
system") within a bipartite quantum system subject to 
Schrodinger evolution. Such effective relaxation to equi- 
librium is predicted to occur under rather general con- 
ditions. Nevertheless deviations are not excluded. Here 
we have been interested in (spatially) selective system- 
environment couplings and in the effect of further parti- 
tioning of the system under consideration. 

For this purpose we have studied a spin-chain consist- 
ing of three interacting spins coupled locally to a quan- 
tum environment. Two different spin-spin-coupling types 
have been examined: a random and a Heisenberg cou- 
pling. For both we have solved the exact Schrodinger 
equation and analyzed the reached equilibrium states. 
We showed that the spin system relaxes into a thermal 
equilibrium state, which is in accordance with the state 
one would expect from quantum thermodynamics. Tem- 
poral fluctuations persist, though, and are a result of the 
still comparatively small environment. 

The spin chain always reaches the predicted state in- 
dependently of the internal coupling type and strength. 
Thus we can conclude, that the internal structure of a 
small system as well as locality of the coupling to a quan- 
tum environment does not effect the relaxation behavior. 
The results also show that the overall state reached is a 
canonical one and thus the spectral-temperature of the 
total spin chain can be identified as the "real" (thermo- 
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